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Abstract 

We calculate the fluctuation of the energy of a system in Tsallis statistics follow- 
ing the finite heat bath canonical ensemble approach. We obtain this fluctuation 
as the second derivative of the logarithm of the partition function plus an ad- 
ditional term. We also find an explicit expression for the relative fluctuation 
as related to the number of degrees of freedom of the bath and the composite 
system. 
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1 Introduction 



Tsallis' statistics, or generalized statistics, represents today a very rich field 
of research. Since its proposal in 1988 by Tsallis B, many developments and 
applications were made in order to give it a physical background. Some examples 
are found in studies of Levy flights (|, |j| , particle physics j|, [| , and turbulence 
[pi . The interpretation relates the parameter q to the fluctuation of some other 
important physical property Hence, many researchers consider it as the 

right framework to describe these phenomena. 

We follow the canonical ensemble with the finite heat bath approach. This 
approach was first proposed by Plastino and Plastino |^| and further confirmed 
by Almeida Q, who gave it a strong mathematical background. In this last 
reference, the author introduced the equation: 



where 0i is given by the Boltzmann principle allowing the calculation of the 
q parameter of Tsallis' statistics, for any system, isolated or not. We adopt a 
convention in the present work: since we are dealing with a canonical ensemble, 
we have three systems, the composed, the studied and the heat bath. Then, any 
quantities pertaining to these systems are labeled by 0, 1 and 2, respectively. 
This is a natural consequence of the structure function formalism |l0| , and its 
connexion to Tsallis' statistical mechanics can be better understood in the fol- 
lowing references, where some issues were treated recently: jyj |l^, 
In this work, another issue is tackled. We propose to generalize the calculations 
of the energy fluctuation within the canonical ensemble. Since the interpretation 
of the parameter q relies on fluctuations, we wish to calculate the fluctuation 
of an important quantity in a canonical ensemble, the system's energy E\ . But 
first, we will revisit the main results in the ordinary picture. 
In the Boltzmann-Gibbs statistics, particularly the canonical ensemble approach, 
we have a general composed system and we wish to study just a small part of it. 
The rest of the composed system is called the heat bath, and it is responsible for 
maintaining the temperature constant. Since the temperature is a constant, so 
is the small part average energy, but the system energy still fluctuates around 
the mean value, and the quantities defined previously can give an estimate of 
these fluctuations. If they are sufficiently small, we can have great confidence 
in the theory results. 

The dispersion of the energy is given, in the canonical ensemble, by the following 
relation, Q: 

. ^ 9 d 2 lnZ 
< A^ 2 >= W (2) 

where Z is the partition function and (3 = -^r anc ^ ^ A-E 2 >=< E 2 > — < 
E > 2 . The relative fluctuation reads, JIo|, |l5|,|l 

<Atf>* 1 (3) 



< E > y/N 
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Therefore, we can conclude that if the number of particles of the system is very 
great, then its energy will be almost all the time very close to its mean value 
< E >, hence confirming the validity of statistical mechanics. 
We will show that, in the generalized picture, the numbers of degrees of freedom 
of the heat bath and of the composite system play a more important role in the 
measure of these fluctuations. 

Another important result will be related to the derivative of lnZ\. We will show 
that the dispersion of energy is still given by equation (||) , but now we have the 
appearance of other term, proportional to < E\ > 2 . 

All these developments will be made in section 2. In section 3, we will illustrate 
these ideas with a very simple situation, two ideal gases in equilibrium. Finally, 
we will address our conclusions. 



2 Calculation of fluctuations 

In Tsallis' picture we have a different energy probability density, which is not 
given by an exponential function, but by a power law in the form: 

P{E X ) = i- [1 - (q 2 - qpEiflte-V firOEx) (4) 

where the parameter (3* is related to the physical temperature by: 

9o-l 1 



q k B T 



and we put a star on it to stress that it is not the Lagrange multiplier associated 
with the entropy optimization problem, as stated in fll7| . Z\ is the partition 
function, and Qi(Ei) is the system's structure function. In general the structure 
function is defined as: 

Qi(Ei) — — l —CiEl l 1 , 
Qi-l 

where the label i represents the i-system 

As was already stated, we do not use g-average values we choose the ordi- 
nary one, defined by the following relation: 

< f(Ei) >= J f(E 1 )P(E 1 )dE 1 (5) 

where this integral, and all others to come, unless otherwise noted, must be 
calculated with the condition on the energy that avoids the distribution P(E\) 
to have negative values, the Tsallis' cut-off condition. This condition has a 
simple physical interpretation in this approach: it means that the system's 
energy cannot be greater than the energy of the whole composite system, that 
is£i<£ , !■ 

Let us calculate the energy mean-square fluctuation of the energy. First of all, 
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we know that the temperature, ksT, is proportional to the energy per degree 
of freedom of a system [|ll] : 

k T = E ° = <El> = <E i> 
B «b/(flb-l) 9i/(9i-l) 92/(92 -1) lJ 



This is the well known equipartition theorem 10 . 

Hence, the energy second moment is calculated as follows: using (||), the relation 
for [3* and the system structure function: 

i i L0 

52-1 p — +\ 



<Ef>= J - { 1 1 



j{i-%y^Er*dE X 

These integrals yield beta functions: 
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<EI>=EI^Z1 (7) 

( gl 12 \ 

\qi-l-> 92-1 J 
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With the knowledge of this quantity, we are able to calculate the energy disper- 
sion. By using the properties of the beta functions and equation (|^), we obtain 
for < (A£i) 2 >: 

. f AP . go/(gO - 1) gl/(gl - 1) + 1 , p . 2 . p .2 

< (Ai/ij >- — r- — — — — - < E x > -<E X > 

- 1 9o/(<Zo - 1) + 1 



Using the additivity of the degrees of freedom: 

9o Qi 92 



9o-l 92-1 
we obtain the relative fluctuation of E±: 



< (AgQ 2 > L 2 
<Ei> 



92/(92 - 1) l 

90/(90 - l) + 1 91/(91 - l) 



(8) 



(9) 



Comparing to equation (^]) we see that it still preserves the inverse square root 
dependence in the number of components of the studied system, since the quan- 
tity 91/(91 — 1) is proportional to it. But, now, this constant of proportionality 
is not trivial. We have a dependence on the degrees of freedom of the whole 
system and the bath, represented by the 90 and 92. Although this happens, we 
still arrive at relation (j^) in the limit 92 — > 1 since, by equation (|§|), when the 
above mentioned limit is taken, the limit 90 — > 1 is implied. 
There is another way of calculating this relative fluctuation, which is to relate 
it to the system's partition function Z\. This is the standard procedure in 
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Boltzmann-Gibbs (BG) statistics. We start with the definition of the energy 
second moment: 



We know that: 
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?2-l ?o 1 d f ( q - 1 E x \ 



< E >= ™ - ~ Ml - E l n 1 {E l )dE 

q 2 9o-l Zi d(l/k B T) J \ q Q k B T 1 

This expression can be rewritten as: 

<EI>=-P- 1 qa 1 9 



n 



q 2 go -l Z l d{l/k B T) 



q - 1 Ei \ / go - 1 Si 



52 — 1 



9o k_bT/ V 90 k B T J 
Using equations (|^) and (^|), we arrive at the dispersion of the energy: 

,/ap^2._ d 2 lnZ 1 9i/(9i-l) + l p 2 (m) 

<{AEl) >~ WIW ~ 9!/(9! - D[9o/(9o - 1) + 1] K El > (10) 

We notice the difference between this result and the usual BG one, the last 
term on the right hand side of equation ( |Io| ) . When we take the infinite bath 
limit, 92 — > 1, the total number of degrees of freedom goes to infinity as well, 
by relation (@) qo —> 1, and this last term vanishes. 

It is a simple algebra task to show that this equation leads to relation (^) for 
the system's relative fluctuation, assuming the equipartition of energy. 
It was shown |]l3| that the usual mean values have the same relations with the 
partition function within this approach and the BG one. All of them are related 
in the same fashion to lnZ\ and its derivatives. Here, the dispersion of energy is 
still related to the second derivative of InZi, but another term comes up, which 
is a function of the degrees of freedom of the system and of the composite one. 
Why we have such a difference? The answer lies in the fact that this approach 
and the BG one have identical energy first moments, equipartition of energy, 
but this fact does not repeat in higher moments. This can be shown by a direct 
calculation of the nth moment of the energy within both approaches. We will 
start with the Tsallis formalism. 
The nth moment of the energy is given by: 



/ (l-f ) n - X E^ 1 {E x )dE 1 
< E[ l >= — i ° J ; (11) 

/(l-f)* 2 - 1 ^)^ 
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It readily gives: 



D 
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<E?>=EZ^^ (12) 

B | -21 22-) 

D \qi~l 1 92-1 ) 

Using the properties of the beta function, we write this expression as: 

, p n >= (gi/(gi - 1) + n - 1) ... (91/(91 - 1) + 1) gi/fa - 1) 
1 (W(<?o - 1) + n - 1) ... (g /(9o - 1) + 1) 9o/(9o - 1) 
In the BG framework, we define the mean value of a function as: 

fe-f> E fSl(E)dE 
1 fe-P E n{E)dE 

Here, we do not need to worry with the cut-off condition, since the whole energy 
is infinite, according to the BG approach. Then, the nth moment of the energy 
is given by the next equation: 



< E n >-- 



1 F U^T 



Once again, using the properties of the gamma function and the value of /3, we 
have: 

< E n >= (91/(91 - 1) + n - 1) ... ( qi /( qi - 1) + 1) 91/(91 - \){k B T) n (13) 



Indeed, equations ([LI]) and (|13| ) have equal energy first moments. In fact, for 
n = 1 in ([Tl|), we have: 

^ I? v. gl ^0 gl , rp 

< E l >= -— — = -k B T 

91-190/(90-1) 9i-l 

with the help of equation (0). Also, for n — 1, equation (13) yields: 

< Si >= -^-ksT 

9i-l 

We note a close resemblance with expression (^2j). In fact, the nth power of 
fcflT is the approximate expression of the product: 

E n 

/jr lim 

1 B ' 9o™ (90/(90 - 1) + n - 1) ... (90/(90 - 1) + 1) 9o/(9o - 1) 

This happens because when we take this limit, the quotient 90 ~ becomes much 
larger than the greatest value of n in the product, provided we choose n not to 
be of the magnitude of the quotient. 

Now if we allow the number of degrees of freedom of the system to be also great, 
much larger than the number n — 1, we end up with the final moment in both 
approaches: 

< E[ l >=< E t >" , 

recovering the equivalence between the canonical and the microcanonical en- 
sembles. 
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3 Application 



To illustrate our results, we take a system which is composed of two ideal gases. 
This composed system has ^ degrees of freedom, and each component gas 
has -^j and degrees of freedom. Before we proceed, a word must be said 
about the treatment of this system in Tsallis' framework. Within our approach, 
any system which is in thermal equilibrium with a finite heat bath, i. e., a bath 
where q 2 7^ 1 , has a probability density given by equation (S) . Hence we do not 
need to have non-trivial effects, such as long range interactions or fractality for 
a system to exibith power law behavior. There is a common belief that these 
last effects are responsible for such framework to arise. 

Our approach relies on another fact, the finiteness of the heat bath. Therefore, 
the main point in our framework is that any system in thermal equilibrium with 
a heat bath which have a finite number of degrees of freedom, q 2 ^ 1, has an 
energy density given by Tsallis power law function. If we take two finite ideal 
gases and make them interact and achieve the equilibrium situation, then we 
can use Tsallis statistics to study this problem. 

We know, from the geometry of the phase space and equation ([!]) , that the 
following relations are valid for an ideal gas: 

« ^ (14) 
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where JV» is the number of particles in the z-th gas and D is the spatial dimension. 
Using equation (^), we have: 

<(AE 1 f>i ( N 2 \ h 1 



<E 1 > ^JVo + iy VnI 

which is still proportional to the inverse square root of the number of particles 
in the studied gas. It is important to strees two properties of this relative 
fluctuation. 

The first one regards the infinite bath limit, N 2 — ► 00. When we take this 

limit, it is implied that the total number of particles oes to infinity as well, 

1 

iVo — * 00. Then, the proportionality constant will be reduced to (-^) 2 , which 
is the constant obtained in the BG context. 

The second regards the magnitude of the number of particles. If Ni and N 2 are 
of the same order of magnitude, this relative fluctuation will be proportional to 
the inverse of the total number of particles. Therefore, it will only be small if 
this total number of particles is large. For example, if we have the numbers: 
Ni = 2, N 2 = 4 and D = 2, the fluctuation of energy will be given by: 



< (A£x) 2 >i (2 



<E!> \7 



0.534 



This means that when we decrease the number of degrees of freedom of the 
thermal bath the energy flutuaction increases. As another exemple if we take a 
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system with N x = 2000, N 2 = 4000 and D = 2 we obtain: 

v J;' 0.018, 

< Ei > 

and now the energy fluctuation is small. 

4 Conclusions 

We showed that the relative fluctuation of the energy of a system in thermal 
equilibrium with a finite heat bath is given by a product of the number of de- 
grees of freedom of the system, the heat bath, and the whole composite system. 
We also showed that the dispersion of the energy still have a dependence on the 
second derivative of lnZ\, but now a new term arises. This term is proportional 
to the square of the mean energy and vanishes properly in the infinite bath limit 
92 -» 1- 

These results show that, within this framework, the role played by the heat 
bath is more explicit than in the BG one. In this picture, the bath has the 
function of keeping the temperature constant and no other quantities related to 
it enter into the equations of the system with which it is in contact. Here, the 
number of degrees of freedom of the bath must be known in order to calculate 
the fluctuations of the energy. This fluctuation cannot be considered small for 
any system, only for large ones when this is an obvious fact. 
We treated an example here, two finite ideal gases in themal equilibrium. The 
relative fluctuation of the system energy depends on the number of particles of 
the bath, iV 2 , and of the composed system, N n . This fluctuation reduces to the 
result in BG statistics in the infite bath limit, N 2 — > 00. Also, if the number 
of particles of both gases are of the same magnitude, the relative fluctuation is 
only small if the total number of particles is large. 
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